Superconductivity in a Hubbard-Frohlich Model and in cuprates. 
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Using the variational Monte-Carlo method we find that a relatively weak long-range electron- 
phonon interaction induces a d-wave superconducting state of doped Mott-Hubbard insulators 
and/or strongly-correlated metals with a condensation energy significantly larger than can be ob- 
tained with Coulomb repulsion only. Moreover, the superconductivity is shown to exist for infinite 
on-site Coulomb repulsion, removing the requirement for additional mechanisms such as spin fluctu- 
ations to mediate d-wave superconductivity. We argue that the superconducting state is robust with 
respect to a more intricate choice of the trial function and that the true origin of high-temperature 
superconductivity lies in a proper combination of strong electron-electron correlations with poorly 
screened Frohlich electron-phonon interaction. 
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It is now over 20 years since the discovery of the first 
high temperature superconductor— and yet, despite in- 
tensive effort, the origin of the superconductivity remains 
fundamentally unknown, with no widely accepted theory. 
The absence of consensus on the physics of the cuprates 
and the recent discovery of iron-based compounds with 
high transition temperatures has re-emphasized the im- 
portance of understanding the origins of superconductiv- 
ity in quasi-two-dimensional materials 2 -. The canonical 
BCS-Migdal-Eliashberg theory could hardly account for 
well documented non-Fermi-liquid properties of cuprate 
superconductors. Moreover, calculations based on the lo- 
cal density approximation (LDA) often predict negligible 
electron-phonon interaction (EPI) insufficient to explain 
a kink in the quasiparticle energy dispersion observed by 
ARPES 3 . Hence, it is not surprising that a large number 
of researchers held the view that the repulsive Hubbard 
model would have the essential physics to account for 
the superconducting and non-Fermi-liquid normal states 
of cuprates. The idea behind this, originally proposed 
by Andersons, is that mobile hole pairs are created via 
a strong on-site repulsion, U. Results by Paramekanti 
et al. — and Yamaji et al4, using a variational Monte 
Carlo (VMC) simulation with a (projected) BCS-type 
trial wave function, appeared to back this up. 

However recent studies by Aimi and Imada£, using a 
sign-problem-free Gaussian-Basis Monte Carlo (GBMC) 
algorithm 8 have shown that these variational methods, 
as well as some other approximations, overestimated the 
normal state energy and therefore overestimated the con- 
densation energy by at least an order of magnitude, 
so that the Hubbard model does not account for high- 
temperature superconductivity cither. This remarkable 
result is in line with earlier numerical studies using 
the auxiliary-field quantum (AFQMC)2 and constrained- 
path (CPMC)i£ Monte-Carlo methods, none of which 
found superconductivity in the Hubbard model. 

On the other hand compelling experimental evi- 
dence for a strong EPI has arrived from isotope 
effects^, high resolution angle resolved photoemission 
spectroscopies (ARPES) 12 , a number of optica l 13 i 14 ' 15 , 



neutron-scatterin; 
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; 16 i 17 and some other spectroscopies 
of cuprates^. Here we show that a long-range EPI 
combined with the Hubbard U provides sizable su- 
perconducting order in doped Mott-Hubbard insulators 
and/or strongly-correlated metals. Previous numerics 
have shown that d-wave order could exist in basic models 
of electron-phonon interactions as a consequence of sec- 
ond order effects (which look similar to those in repulsive 
models), with d- selected over s-wave order by the Hub- 
bard repulsion, which acts as a kind of filter—. Other 
work determining the properties of strongly-coupled po- 
larons and bipolarons have shown that the long-range dis- 
crete Frohlich EPI in the presence of a strong Coulomb 
repulsion, does not lead to the enormous enhancement of 
the carrier effective mass characteristic of the Hubbard- 
Holstein model (HHM)22i2L22i2i2i^Si^. It is our 
thinking that the special properties of the Frohlich type 
interaction are also likely to be interesting in the weak 
electron-phonon coupling regime with intermediate car- 
rier density. Specifically, the purpose of this study is to 
determine if d-wave superconductivity is enhanced by a 
combination of the lighter carrier masses and propensity 
to intersite pairing found in models with long-range in- 
teractions and significant, repulsive Hubbard U. 

Our Hubbard-Frohlich model (HFM) contains the 
usual electron hopping, t(n), and electron-electron on- 
site repulsive correlations, U. In addition there is a term 
to describe lattice vibrations, with frequency u> and mass 
M, and a term to describe the long-range interaction of 
the electrons with ion displacements, 

H= f ( n ~ n ')cl v c nCT + U^2 n nT n ni + 

n,n',cr=T,l n 

E \B + - £ !McLc^ m - (i) 



Here c n , a and c n<T create and annihilate the electron with 
spin a at sites n' and n, respectively, P m = —ifid/dt; m 
is the ion momentum operator at site m, and £ m is the 
ion displacement. The long-range Frohlich EPI is char- 
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acterized by the force which has the formS 2 -, 

... k ( |m — n|\ 

/m(n) ^ [(m-n)V^ + l]3/ 2 eXP (-^^J' (2) 

where m, n are the lattice vectors, a is the lattice con- 
stant, and R sc is the screening radius. That is, the 
screened force is the unscreened force multiplied by an 
exponential screening factor. The Frohlich EPI, Eq.0, 
routinely neglected in the Hubbard U and t-J models of 
cuprate superconductors 29 is of the order of 1 eV as es- 
timated from optical data^. The force function, Eq.([2]), 
describes the EPI with c-axis polarized optical phonons, 
which is poorly screened since the upper limit for the 
out-of-plane plasmon frequency in cuprates^ is well be- 
low the characteristic frequency of optical phonons. 

The carrier mass and the range of the applicability 
of analytical weak and strong-coupling expansion effec- 
tively depend on the EPI radius. In particular the exact 
carrier mass calculated with the Frohlich EPI using the 
continuous-time QMC algorithm 2 ^ was found to be sev- 
eral orders of magnitude smaller than in the Holstein 
model in the relevant region of fuv/t ratio. The mass is 
well fitted by a single exponent remarkably close to those 
obtained using the Lang-Firsov transformation and sub- 
sequent averaging over phonons, for any strength of the 
Frohlich EPI, justifying our use of the Lang-Firsov trans- 
formation in this work. Here we use the transformation 31 
to integrate out phonons (for technical details see for ex- 
ample Ref. 28 ) mapping the electron part of the trans- 
formed Hamiltonian, H, to an extended U-X Hubbard 
model with renormalized hopping integrals, t(n), a di- 
minished on-site repulsion, U, and a long-range effective 
attraction, XW, 

H = ^2 t(n - n')4, ff cw + U h n] n ni 

n^n',cr n 

- XW Hn-n')n n , a h na , (3) 

Here t(n) = t(n) exp[-c/ 2 (n)], g 2 (n) = [E p - 
\W$(n)]/hu, $(n) - «- 2 E m /m(0)/ m (n), U = U - 
2E p , W — zt(a) is the renormalized half-bandwidth, 
z is the coordination number and A = k 2 /2Muj 2 W is 
proportional the conventional BCS electron-phonon cou- 
pling constant. The latter is proportional to the single- 
particle density of states (DOS) and could be larger than 
A due to the van-Hove singularity of DOS and correlation- 
enhanced effective mass of carriers. The polaronic shift, 
E p = (k 2 /2Mw 2 )$(0), of atomic levels is included in 
the chemical potential. In the following only nearest- 
neighbor hops are allowed. As shown by Bonca and 
Trugman 2 ^ the physical properties of (bi)polaronic car- 
riers depend predominantly on the EPI coming from the 
first two sites in Eq.Q, so that we take R sc = oo in our 
simulations. 

When EPI is strong compared with the renormalized 
kinetic energy, A ^> 1, one can apply the 1/A perturba- 
tion expansion reducing the multi-polaron problem to a 



charged Bose-gas of small mobile bipolaron o 20 ' 28 . Inter- 
mediate and weak coupling regimes, A < 1, with large U 
require a variational approach. Here we use a standard 
VMC method as for example in Ref. [H to minimize the 
energy H\^f)/\^f). The difference is that we make an ad- 
ditional measurement of the long-range attraction in the 
extended Hubbard Hamiltonian, Eq.Q. 

A BCS type trial wave-function is used, which has the 
form 

l* T > = p n p g JJ( Uk + vk4 T cL kl )|o), (4) 

k 

where Pjv = &y^.hi,N projects onto a fixed particle num- 
ber state and Pq = IlrJl — (1 — flO^nT'W] P r °j ec ts 
out double occupancy 6 -. Our variational parameters are 
g ^ 1 of the Gutzwiller double occupancy pro- 
jection operator, Pq, and the chemical potential, p, 
that enters through the kinetic energy, £k = t(cosk x + 
cosk y ) - p, via tt 2 . = (1 + Ck/VH + A 2 .)/2 and v 2 . = 

(1 — £k/ \/£k + A k )/2. We then vary the parameters to 
minimize the energy for values of the electron density p, 
U, superconducting order parameter Ak and A keeping 
t = 1. We place 42 spin- up and 42 spin-down electrons 
on a 10 x 10 2D square lattice to simulate optimal doping, 
keeping the electron density fixed at p = 0.84. Two cases 
are investigated: a) an on-site repulsion U = 8 is used 
for comparison with results of Ref. @, and b) an infinite 
U is used to see if the attraction induced by EPI, alone, 
is enough for superconductivity. 

We examine the d-wave Ak = A(cosfc :E — cosk y ), 
extended s-wave Ak = A(cosfc a; + cosk y ) and s-wave 
Ak = A order parameters. Following previous work 6 
periodic boundary conditions are used in one direction 
and anti-periodic boundary conditions are used in the 
other. The value of A is varied to study the effect of in- 
creasing EPI and we find the value of A at which the 
energy minimum occurs. We also measure the energy of 
clustered states to ensure our results are stable against 
the formation of immobile clusters. In our results all en- 
ergies quoted are per electron; note that Yamaji et al£ 
use energy per site. 

In the following, all energies are quoted in units of t. 
For the A = and U = 8 case shown in FigQ] (top), we 
recover Yamaji's result 6 with a minimum at A w 0.08 
and energy per electron of (—0.8780 ± 0.0004), the nor- 
mal state energy is (—0.8763 ± 0.0004). It can clearly be 
seen that the effect of increasing EPI is to increase the 
depth of the minimum and therefore the stability of the 
superconducting state. Our minimum energy per elec- 
tron is (-4.4008 ± 0.0004) at A = 0.16 and A = 0.075. 
The maximum condensation energy gain l-E'(A) — i£(0)| 
with A = is (0.0017 ±0.0006) and (0.0049 ±0.0006) per 
electron with A = 0.075. 

A surprising result is found for infinite Coulomb repul- 
sion, where mapping to the t-J model would result in 
J| = At 2 /U = (so no spin fluctuations are present). 
FigfT] (bottom) shows an increasing condensation energy 
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FIG. 1: Condensation energy per electron (in units of t) versus 
the amplitude of the superconducting d-wave order-parameter 
for U — 8, top, and U = oo, bottom, with different EPI 
coupling, A. 



gain with A for U — oo: the maximum condensation en- 
ergy gain is (0.00050 ±0.00007) per electron at A = 0.06. 
While this energy gain is not large enough to guarantee 
stability against other trial wave functions, it does sug- 
gest a d-wave superconducting state exists for U = oo 
with strong enough A, that is d-wave order can be in- 
duced without spin-fluctuations (as expected, the A = 
curve shows no superconducting order). 

We examined the energies of various cluster forma- 
tions. The lowest clustering energy for A = 0.075 is 
—3.6824 per electron, so that the system is stable against 
the formation of clusters. It is found that there is no s- 
wave or extended s-wave superconducting state for both 
U = 8 and U = oo with A up to 0.075 see FigHJ 

Accurate QMC studies have proven that the Hubbard 
U and t-J models are unlikely to account for the high 
critical temperatures of cuprate superconductors 7 . On 
the other hand ab-initio LDA calculations of EPI have 
led to a conclusion that phonons fail to explain high- 
temperature superconductivity either 3 . It is well known 



FIG. 2: Condensation energy per electron (in units of t) versus 
the amplitude of the superconducting s-wave order-parameter 
(top) for U — oo and U — 8, showing no s-wave ground 
state. The condensation energy versus the amplitude of the 
extended s-wave order-parameter (bottom) for U = oo and 
U — 8, showing no extended s-wave ground state. 



that LDA badly underestimates the role of the Coulomb 
correlations, incorrectly predicting that the parent com- 
pounds of cuprate superconductors (such as La2Cu04) 
are metallic rather than insulating. Also, the anisotropy 
of the electron response functions determined with LDA 
is much smaller than the experimentally observed value 
in these layered materials. It is not surprising at all that 
EPI turns out to be very weak in this "metallic" picture 
due to electron screening of the long-range electron-ion 
interaction. Moreover, it has been noted that the inclu- 
sion of Hubbard U via the LDA+U scheme significantly 
enhances the EPI strength 32 since the system becomes a 
doped Mott insulator with poor screening in at least the 
c-axis direction as anticipated in Refs. HH andl20l. 

These theoretical results, as well as many 
experiment o 11 ! 12 ' 13 ' 14 ^ 5 ' 16 ! 17 , tell us that any realis- 
tic theory of high-temperature cuprate superconductors 
should include strong electron-electron correlations and 
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the finite-range EPI on equal footing, as previously 
discussed by Alexandrov and Mott 33 . When an ex- 
tended Frohlich EPI is considered, the picture changes 
drastically compared with the local Holstein interaction, 
where carriers are readily self-trapped on a single lattice 
site even by a moderate electron-phonon coupling 2 ^. It 
is interesting to note the profound effect of long range 
electron-phonon coupling on the magnitude of the su- 
perconducting order, since such long range interactions 
have equally dramatic effects when electron-phonon 
coupling is strong. For example, we found that small 
mobile bipolarons could exist when electrons interact via 
a strong finite-range Frohlich EPI with high-frequency 
optical phonon o 20 ! 25 ' 26 . 

Here, we showed using the VMC method, that even 
a relatively weak Frohlich EPI is sufficient to induce a 
d-wave superconducting state with substantial condensa- 



tion energy in a doped Mott-Hubbard insulator and/or 
strongly-correlated metals. The calculated condensation 
energy and the difference between clustering and super- 
conducting energies in a wide region of A are larger than 
any difference between VMC and the more exact GBMC 
algorithm of Aimi and Imada*, so that our supercon- 
ducting state is very likely to be robust with respect to a 
more complicated choice of the trial function. As a result 
we believe that the true origin of high-temperature su- 
perconductivity lies in the proper combination of strong 
electron-electron correlations with long range electron- 
phonon interactions. 

We would like to thank Pavel Kornilovitch for long- 
standing collaboration and helpful discussions. We ac- 
knowledge support of this work by EPSRC (UK) (grant 
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